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1 Introduction 



The set of orthogonal polynomials P n {x) = x n + ■ • ■ , with respect to an exponential weight 



/oo 
P n {x)P m {x)e v{c > x) dx = h n 5 nm , V(c,x) :=J^c fc 
-°° k>i 



x h , 



is an essential ingredient of the methods pQ -[2] for studying the large-n limit of the Hermitian 
matrix model 

Z n = JdMexp (TrV(c,M)), (1) 
One of the main tools used in these methods is the pair of equations 

zP n (z) =ZP n (z), d z P n (z)=MP n (z), n>0, (2) 
where (Z, Ai) is a pair of Lax-Orlov operators of the form 

Z = K + u n + v n K\ M = ~Y,kc k (Z k - 1 ) + . (3) 

k>l 

Here A is the shift matrix acting in the linear space of sequences, A* is its transposed matrix and 
( )+ denotes the lower part (below the main diagonal) of semi-infinite matrices. 

The first equation in ([2]) represents the standard three-term relation for orthogonal polynomials. 
Both equations are referred to as the string equations in the matrix models of 2D quantum gravity 
m and provide the starting point of several techniques to characterize the large-n limit of (pQ). A 
deeper mathematical insight of these methods was achieved after the introduction by Fokas, Its 
and Kitaev [3] of a matrix valued Riemann-Hilbert (RH) problem which characterizes orthogonal 
polynomials on the real line, and the formulation by Deift and Zhou [2],[1] of steepest descent 
methods for studying asymptotic properties of RH problems . 

The RH problem of Fokas-Its-Kitaev was generalized by Van Assche, Geronimo and Kuijlaars 
[5] to characterize multiple orthogonal polynomials. Moreover, it was found [6] -[TO] that these 
families of polynomials are closely connected to important statistical models such as Gaussian 
ensembles with external sources and one-dimensional non-intersecting Brownian motions. 

In this paper we generalize the string equations ([2]) to multiple orthogonal polynomials of types 
I and II, and show how these equations can be applied to analyze the large-n limit of multiple 
orthogonal polynomials and their associated statistical models. Section 2 introduces the basic 
strategy of our approach to derive string equations, which is inspired by standard methods used 
in the theory of multi-component integrable systems [llj-|15|. As it was proved in [H] the multiple 
orthogonal polynomials of types I and II are elements of the first row of the fundamental solution 
/ of the corresponding RH problem. Then, in Sections 3 and 4 we formulate systems of string 
equations for the elements of the first row of the fundamental solution /. In both cases the function 
/ depends on a set of discrete variables 



(si, S2, . . . , s q ) e Z q , where 



Si > for type I polynomials 
Si < for type II polynomials. 
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Therefore, special care is required to determine the form of the string equations on the boundary 
of the domain of the discrete variables. Thus, we obtain closed-form expressions, free of boundary 
terms, for the string equations satisfied by these types of multiple orthogonal polynomials. These 
string equations are associated to pairs (Z{,A4i) of Lax-Orlov operators. In particular those 
involving the Lax operators Z{ lead to the well-known recurrence relations for multiple orthogonal 
polynomials [5]. 

We take advantage of an interesting observation due to Takasaki and Takebe [16] who showed 
that the dispersionless limit of a row of a matrix-valued KP wave function is a solution of the 
zero genus Whitham hierarchy [11]. This is an additional incentive for using Lax-Orlov operators 
|12]-[15] in order to characterize the large-n limit in terms of quasiclassical (dispersionless limit) 
expansions. Thus, in Section 5 we show how the leading term of the expansion of the first row 
of / is determined by a system of dispersionless string equations for q + 1 Lax-Orlov functions 
(z a ,m a ) in the phase space of the Whitham hierarchy. The unknowns of this system reduce to a 
set of q pairs of functions (uk,Vk), which are determined by means of a system of hodograph type 
equations. Finally, Section 6 is devoted to illustrate the applications of our results to models of 
random matrix ensembles and non- intersecting Brownian motions. 

The present work deals with multiple orthogonal polynomials of types I and II only, but the 
same considerations apply to the study of multiple orthogonal polynomials of mixed type [T?] • On 
the other hand, we concentrate on the description of the leading terms of the asymptotic solutions 
in the dispersionless limit. However, as it was showed in [18j-|19j for the case of the Toda hierarchy 
and the Hermitian matrix model, the scheme used in the present paper can be further elaborated 
for determining the general terms of these expansions, as well as their critical points and their 
corresponding double scaling limit regularizations. 

2 Riemann-Hilbert problems 

In this work we will consider (q + 1) X (q + 1) matrix valued functions. Unless otherwise stated 
Greek a, (3, . . . and Latin i, j, . . . suffixes will label indices of the sets {0,1,..., q} and {1,2, ... ,q}, 
respectively. We will denote by E a p the matrices {E a p) a >p> = 5 aa > 6/3/31 of the canonical basis and, 
in particular, its diagonal members will be denoted by E a := E aa . Some useful relations which 
will be frequently used in the subsequent discussion are 

E a f3 E^x = <5/3 7 E a x; E a aEp = a a p E a p, V matrix a. 

We will also denote by V(c, z) the scalar function 

V(c,z) :=Y J C n z n , c=( Cl ,c 2 ,...)GC 00 , (4) 

n>l 

and will assume that only a finite number of the coefficients c n are different from zero. 

Given a matrix function g = g(z) (z E R) such that det g(z) = 1, we will consider the RH 
problem 

m_(z) g{z) = m + (z), (5) 

where m(z) is a sectionally holomorphic function and m±(z) := lim e ^o+ m ( z ± ie). We are 
interested in solutions / = f(s, z) of © depending on q discrete variables s = (si, . . . , s q ) E 
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such that 



f(s,z) = [l + O[-))f (s,z), z^oo, 



(6) 



where 

q q 
f (s,z) :=^2z Sa E a , (s := Sl ). 

a=0 i=l 

The set of points s £ Z 5 for which ([5]) admits a solution f(s,z) satisfying ([6]) will be denoted by 
r. The solution f(s,z), (s G T) is unique and will be referred to as the fundamental solution of 
the RH problem ([5]). 

We will apply ([5]) and ([6]) to derive certain difference-differential equations for /. These equa- 
tions contain two basic ingredients: the coefficients of the asymptotic expansion of f(s, z) as z — > oo 



f(s,z) = (l + Y t 



a n (s) 



fo(s,z), 



n>l 



and the q pairs of shift operators Tj, T* acting on functions h(s) (s G T) defined as 
h(s — ej) if s — e-i G T 



(Tj /t)(s) 



x o if s - ei i r 

where are the elements of the canonical basis of C 9 . 
We will often consider series of the form 



h(s + ej) if s + ej G r 
[0 if s + ei i T, 



n=l 



n=l 



and will denote 



ra=l 



n=l 



(7) 



(8) 



The RH problem © admits the following symmetries. 

Proposition 1. 1. Ifh(s,z) (s G T) is an entire function of z, then h(s, z) f(s, z) satisfies © 
for all s G r. 

2. The functions (T, f)(s, z) and (T* f)(s, z) satisfy © for all s G T. 

3. If g(z) is an entire function, then for any entire function (j)(z) verifying 

g- 1 (f)g = <\>-g~ x d z g, 



the covariant derivative 
satisfies ([5]) for all s G T. 



D z f:=d z f-f< 



(9) 
(10) 
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Our strategy to obtain difference-differential equations for / is based on applying the next 
simple statement to the symmetries of ([5]). 

Proposition 2. Let f(s,z) be a a solution of © defined for s in a certain subset Tq C T. If 
f(s, z) f(s, z)~ x — P(s, z) — > as z — > oo ; where P(s, z) is a polynomial in z, then 

f(s,z)=P(s,z)f(s,z). 

Proof. Since det g(z) = 1 it follows from ([5]) and © that det f(s, z) = 1 so that the inverse matrix 
f(s, z)^ 1 is analytic for z £ C — R and satisfies the jump condition 



g(z)- 1 f-(s,z)- 1 =f + (8,z)- 1 , z£R. 
As a consequence / / _1 is an entire function of z and the statements follow at once. 



□ 



3 Multiple orthogonal polynomials of type I 

Given q exponential weights Wi on the real line 

Wi {x) := e~ V ^ x \ Cl = (c a , Ci2, ■ ■ ■) € C°°, 
and n = (m, . . . , n q ) E N 9 with |n| > 1, the type I orthogonal polynomials 

A(n, x) = (Ai(n, x),..., A q (n, x)) 
are determined by means of the following conditions: 

i) If rij > 1 then the polynomial Aj(n, x) has degree rij — 1. If rij = then Aj(n, z) = 

ii) The following orthogonality relations are satisfied 



dx 



l = 0,l,...,|n|-2, 

1 Z = Inl - 1. 



We assume that all the multi-indices n are strongly normal [T7] so that A/(n, z) are unique. 
The RH problem which characterizes these polynomials j5] is determined by 



9(z) 



( 1 ... \ 

-wi (ar) 1 ... 



\ -w q {z) ... 1 / 
The corresponding fundamental solution /(s, z) exists on the domain 

r ; = {seP : > 0, Vi = 1,... 



(11) 



(12) 
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For s ^ it is given by 



f(s,z) 



R(s,z) 



( R(s,z) A(s,z) \ 

d^R(s + ei, z) d^ 1 A(s + e\,z) 

\ d q l R(s + e q , z) d q 1 A(s + e q ,z) J 



dx Y?j=i A 3( a > x ) w j( x ) 



2ni z — x 

where dj is the leading coefficient of Aj(s + ej, z). Furthermore, for s = 



f(0,z) 



( 1 

Ri{z) 1 
R 2 (z) 1 



V Rq{z) 

Because of the form of Ti we have that 



\ 



1 ) 



Rj(z) 



dx Wj{x) 
2iri z — x 



h(s — if Si < 1 



(T*h)(s) = h(s + ei ), (Tih)(s):-- 



(13) 



(14) 



if 5, = 0, 

for functions h(s) (s 6 Tj). It is clear that 

T*T l = I, TiT* = {l-6 Si>0 )I, 

where I stands for the identity operator. Sometimes it is helpful to think of the functions h(s) 
as column vectors (h\ Si= o, h\ Si= \, h\ Si =2, ■ ■ -) T - Thus, in this representation, Tj, T* become the 
infinite-dimensional matrices 

/ 1 \ 

10 
1 



V 








... 


" ^ 


1 













1 







V '■■ 









3.1 The first system of string equations 

From the asymptotic expansion ([7]) we have that as z — > oo 



zE + a(s - E - E a(s) + I - E - E { + 0[ -j , Vs G T T + e u 
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where we are denoting 

a(s) := ai(s). 
Hence by applying Proposition 2 it follows that 



{Ti f)(s, z) = (zE + a(s - a) E - E a(s) + I - E - E t ) f(s, z), seTj + e, 
which implies 

(T l Eof)(s,z) = ((z-u i (s))E -Y,a 0j {s)E 0j ^f(s,z), VseT I + e l , (15) 

j 

where 

Ui(s) := a 00 (s) - a 00 (s - ei). 

Similarly one finds 

(T*E f)(s,z) = a 0j (s + e j )E 0j f(s,z), Vs 6 T/. (16) 
Note that as det f(s, z) = 1 for all (s, z) G Yj x C then , as a consequence of (fl6]l we deduce that 

a oj (s + ej)^0, Vs G T/. 

If we now define 

a j(s + ej) 

then from (I15p it follows that 

Proposition 3. T/ie function f satisfies the equations 

z (E f)(s, z) = (Ti + Ui (s) + Vj(s) T/) (Eo f)(s, z), (18) 

j 

for all s GTi + ei and i = 1, . . . , q. 

As a consequence we get the following system of string equations 
Theorem 1. The multiple orthogonal polynomials of type I verify 

zA(n,z) = (Ti + Ui (n) + Y,Vj(n)T*) A(n,z), (19) 



for all n €Tf + ei and i = 1, . . . , q. 

For q = 1 Eq. (|19[) reduces to the classical three-term recurrence relation for systems of orthog- 
onal polynomials on the real line. 

On the other hand Eq, ()19() implies 

A(n-ei,z)-A(n-ej,z) = (a 00 (n- ei) - a 00 (n- e,-)) A(n, z), Vs £ T I + e i + e j , i ^ j. (20) 

The relations (|19p and (|2Up lead to a recursive method to construct the multiple orthogonal poly- 
nomials of type I. Indeed, it is clear that for n = n^ej we have 

A( ni ei,z) =Af(n h z) ei = (0, . . . ,0, (m, z),0, . . . ,0), 

where A^\ni,z) are the orthogonal polynomials with respect to the weight Wi(x). Then starting 

from Af\rii,z) and using (|20p we can generate the multiple orthogonal polynomials of type I for 
higher q. 
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Example 

Let us denote by Ij n the moments with respect to the weight Wj 



>./.» ■• I ^ n ™^\ (21) 



We have that 

A fl \ 1 A (O \ ^'° Z ~ Ix ' 1 



h,0 ' ' ^1,0-^1,2 — 

The recurrence relation (|19p for q = 1 is 



A x {n + l,z) = Q0l(n + X) [( z + aoo (n - 1) - a Q0 {n))A x {n, z) -A x {n- 1, z)], Vn > 2, (22) 
a i(n) 



where according to (fT3j) 

/* dx 

aoo(n) = / — x n A(n,x)wi{x), (23) 
Moreover, the normalization condition gives us 

'"" ' " ' f ^ n [{x + a 00 (n - 1) - a 00 {n))A{n, x) - A(n - 1, x)] Wl (x). (24) 



aoi(" + 1) 7 27ri 

The system (I22p -(|23p- (|24p allows us to construct the polynomials A(n,z) for n > 3. For example 
one gets 

A , g ^ = (A 2 1 ~ V^g) ^ ~ A,l A,3 + (/l,0/l,3 ~ + A 2 2 

^1,2 ~~ (2-^1,1-^1,3 + h,oh,4)h,2 + hfll\fl + ^1,1^1,4 
If we write (I20p in the form 

A(n — es,z) — A(n — e,-, z) „ . . 

aoo(»-ei)-aoo(n-ej) 

and take into account that 

r dx q 

aoo(n) = / —x\ n \y^A k (n,x)w k (x). (26) 

we can construct all the multiple orthogonal polynomials of type I. Thus, for example for q = 2 we 
obtain 

A(1,M) = ^ (/ 2 ,o,-/i,o), 

A(2, l,z) = — (/i, 2 /2,o - /i,i/2,i + z(h,oh,i - h,ih,o) , - A.o-^o) , 



where 



C 2 

Cl : = A, 1-^2,0 — ^1,0-^2,1) 



C*2 : — ^2,2-^1,1 — ^1,3-^2,0-^1,1 — ^1,2-^2,1-^1,1 + A,2-^2,0 + ^1,0-^1,3-^2,1 — ^1,0-^1,2-^2,2 • 
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3.2 Lax operators 



The functions foi(s, z) = Ai(s, z) can be written as series expansions of the form 



foi(s,z) 



z z z 



where 



OLinis) = 0, Vn > Si + 1. 



On the other hand, it is easy to see that 
Hence, from (p7|) and (f28j) it is clear that 



n-l 



k=0 



—2 ^ -z s ' = a; in +i(s + e J )7f z , Vn>l, 

z 

so that we may write 

/oiO + ei,2f) = z), £i(s,z) := z Si , seTj, 

where the symbols Gj are dressing operators defined by the expansions 

Gi = a in (s + ej) T™ _1 , a in (s + e») := (a n )oi(s + e i), 

n>l 

or, equivalently, by the triangular matrices 



G; 



/ Goo 










" ^ 


Gio 


G11 










G20 


G21 


G22 








G n rn — C^i,?!— m+1 "I - ^i) 



The inverse operators can be written as 



aii(s + ei) a i(s + ej 



We define the Lax operators .Zj by 



Zi := GiT* GJ 1 . 



It follows at once that they can be expanded as 

Z i = li {s)T* + Y< r 1™{s)'lT, 



n>0 



where 



7i (a) = + ei ) (2? #i)00 = a f" x + 9 ei) , =Vi( 8 + ei). 

aa{s + 2ej) 



Proposition 4. The functions /oj satisfy the equations 

z f oi (s + ei,z) = (Zif oi )(s + ei,z), Vs G I 1 /. 
Proof. From the definition of Gi we have 

zf oi (s + e i ,z) = G i (zC i ) = (G i T*)(C i ) = (Z i f 0i )(s + e i ,z). 



3.3 The second system of string equations 

Let us consider diagonal solutions 

${z) = Y J Mz)E a 

a 

of the condition ([9]) corresponding to the function g(z) of (jlip . They are characterized by 

<9 2 - 0o Wi + <fo = 0, i = l,...,q. 
In this way, by setting 4>q = we get 

The corresponding covariant derivative is 

D z f:=d z f-J2v'{c i ,z)fE i . 

i 

Hence we have 

D z ( E f) = d z f 00 E + ^2(dJ 0i - V'( Ci , z) f oi ) E oi . 

i 

It is clear that ([33]) implies 

z n f 0i (s + e h z) = (Zf / w )(s + e u z), Vs 6 !>. 
On the other hand, as z — > oo 



(TO n /0a)(*^)=< 



— )z s °, for a = 0, 



, n > 1, 



) zSi , for a = i / j , 



(2T/oi)(*,«) = ^ 



°(^+Tl **> for Si - n ' 



, n > 0. 



0, for Si < n 
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Proposition 5. The function f satisfies the equation 

(D z +H)(E f)(s,z) = 0, Vser 7 + ^ ei , 



where Ti is the operator 



3=1 



(38) 



(39) 



Proof. Given s £ Tj + ^ • ej let us denote 



From (1351) it follows that 



(I> z + W) (£ /) 



2 /oO + X] V '( C 3> Z j)(j,+) /< 



00 



-En 



(40) 



" 



Now from (|37p we have 
g 



d z foo + ^'(^,2^/00 = 0(- 



i=i 



d z f 0i +Y,v^ J ,z J ) { 



O -)z Si . 



W'^Z^-V'^zj) f 0i (s,z) = (y'(ci,Zi)-V'(ci,z)) f oi (s,z) + o(± 
Moreover, from (j36[) it is clear that 

(y'faZi) - V'(ci, z)) /oi(a, z) = (y'faZi) - V'( Ci , z)) / *(s W + Ci, z) 



^nq R (^f- 1 - z"- 1 ) /oiCaW + e i; z) = 



n>l 

Therefore we find 



{D z + W) (E /)(«, z) = O(-) f (s, z), z -> oo. 



The first member / := (D z + TL) (Eq f) of this equation is a solution of ([5]) for all sG T/ + ^ ■ 
and /(s, z) /(s, z) _1 — > as z — > oo. Therefore, the statement of Proposition 2 implies / = . □ 



As a consequence we deduce the following system of string equations 
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Theorem 2. The multiple orthogonal polynomials of type I verify 

q 

d z Ai(n,z) =V'(c i ,Z i )A i (n,z)-^2v'(c j ,Z j ) ijt+) A i (n,z), (41) 

j=i 

for all neT/l and i = 1, . . . , q. 

3.4 Orlov operators 

We define the Orlov operators Mi by 

Mi'-Gi-Si-Ti-Gr 1 . (42) 
They satisfy [Zi, Mi] = I and can be expanded as 

Mi = J2»in(s)T?, (43) 

n>l 

where 

WW = (44) 
Proposition 6. T/ie functions /oi satisfy the equations 

d z foi(s + ei,z) = (Mifoi)(s,z + ei), Vs G T 7 . (45) 
Proof. From the definition of Gj we have 

9 2 / j(s + e*, z) = d {si z~ l ii) = G- Si{Ti &) = G f • a* • T f ■ G'Voi = -Mi /oi(« + c*, z). 

□ 

4 Multiple orthogonal polynomials of type II 

We consider now q exponential weights Wi on the real line 

Wi (x) := e v <**\ c l = (c il ,c l2 ,...)GC°°. 

Note the difference in the sign of the exponents with respect to the weights for multiple orthogonal 
polynomials of type I. Given n = (n\, . . . ,n q ) £ N 9 , the associated type II monic orthogonal 
polynomial P(n, x) = x'™' + • • • is determined by the conditions 



/ P(n, x) Wi(x) x 3 dx = 0, j = 0, . . . , rtj — 1. 
Jr 



We assume that all the multi-indices n are strongly normal [T7] so that P(n, z) is unique. 
The RH problem for the multiple orthogonal polynomials of type II is determined by [T7] 



9(z) 



( 1 w%(z) w 2 {z) . . . w q (z) \ 
1 ... 



\ ... 1 / 



(46) 
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Its fundamental solution f(s,z) exists on the domain 

r n = {seZi : Si <0, Vi = l,...,g}, 

For S{ < — IVi, it is given by 



(47) 



f(s,z) 



( P(n,z) R(n,z) \ 

d\ P(n — ei, z) d\ R(n — e±, z) 

\ d q P(n- e q ,z) d q R(n-e q ,z) J 



n. 



Rj(n,z) := / 

JR 



(48) 



dx P(n, x) Wj{x) 
2ni x — z 



dx 



P(n — e,j , x) Wj (x) x r 



For the remaining cases, in which one or several Sj vanish, one must insert the following corre- 
sponding row substitutions in (f4~8|) 



(diP(n- ei,z) diR(n- ei,z)) — ► (0 e^. 



(49) 



In particular 



/ 1 R^z) R 2 (z) 
1 



f(0,z) = 

\ 

In view of (|47|) we have that 



RJz) \ 



i 7 



Rj(z) 



dx Wj(x) 
2ni x — z 



{T i h)(s) = h(s-e i ), (T*h)(s):={ 



h(s + if Si < —1 
if Si = 0, 



for functions h(s) (s £ T//). Note also that 

Ti^* = i, i7ri = (i-5 s . )0 )i 

where I stands for the identity operator. If we think of h(s) as a column vector (/i| Sj= o, /i| Sl= _i, /i| Si= _2, • • -) T , 
then Tj, T* are represented by the infinite-dimensional matrices 



/ 1 ... 
10 
1 



3 



V i i '-J 



( ... 
10 
10 



\ 
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4.1 The first system of string equations 

The same analysis as in the subsection 3.1 leads now to the equations 

(T i E o f)(s,z) = Uz-u i (s))E -J2ao j (s)E oj )f(s,z), Vs € T n , (50) 



where 

Ui(s) := a o0) - aoo(« - e*). 

Similarly one finds 

(T*E f)(s, z) = a 0j (s + ej) E 0j f(s, z), Vs G T n - e h (51) 
and taking into account that det f(s, z) = 1 for all (s, z) £ Tjj x C, from (|51 j) we obtain 

a 0j (s)^0, VsgT/j. 

Now we define 



Vj{s) := < 



a j(s + ej) 



(52) 



0, for Sj = 0. 



Notice that the functions Vj(s) (s G T/) for multiple orthogonal polynomials of type I defined in 
(fT7|) also satisfies i>j(s) = for Sj = 0. 
If we now recall that according to (|49p 

/(«, ^) = E 0j , for = 0, 

from (|50p it follows that 

Proposition 7. T/ie function f satisfies the equations 

zE f(s, z) = (r { + Ui (s) + ^ Vj(s) if) {E f)(s, z) + £ a » (53) 

i 3 

for all s £ Tjj and i = 1 . . . q. 

As a consequence we get the string equations 
Theorem 3. The multiple orthogonal polynomials of type II verify 



P(n, z) = (Ti + Ui (-n) + v^-n) T*) P{n, z), (54) 



for all n and i = 1, . . . , q. 

These equation provide a recursive method to construct multiple orthogonal polynomials of 
type II. We may write ([54"]) as 

9 (n) 

P(n + e j ,z)-a 00 (n + e j )P(n,z) = (z-a 00 (n))P(n,z)- V] -P(n-e k ,z), (55) 

k=t^ k >i aok{n ~ &k) 
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where, according to(|48|). we have that 

aoo(n) = coeff[P(n,z),z |n ' l_1 
dx 



aofc(n) 



2m 



P(n, x)x nk w k (x)dx. 



On the other hand, multiplying the equation (|55|) by z n ^>Wj(z), integrating on 
orthogonality condition for P(n + ej, z), we obtain 



(56) 



and using the 



a 00 (n + ej) 



f dx 

/ P(n,x)x nj Wj(x) 

Jr 2m 



(x - a QQ (n))P(n,x) - ^ 



fc=l,nfc>l 



a 0k (n - e k ) 



P(n - e k ,x) 



x nj Wj(x)dx, 



so that 



a 00 (n + ej) 



a 0j (n) 



[x 



a 00 (n))P(n,x)- £ ^ {n - e k ,x 



x nj Wj(x)dx. 

(57) 

The system (I55p - (l56p - (l57p determines the multiple orthogonal polynomials of type II in terms of 
the moments Ij, n . 



Example 

For q = 1 is clear that 



P(0,z) = l, P(l,z)=z-!±±. 



From (f55l) - (f56l) - (f57l) we easily obtain that 



P(2,z) 
P(3,z) 



z 2 + (h,oh,3 - h,ih,2)z h,2 ~ h,ih,3 
I'll — h,oh,2 ^ii - h,oh,2 ' 



3 ( — ^1,5-^1,1 + ^1,3-^1,1 + h, 2-^1,4-^1,1 — -^1,2-^1,3 — h,oh,3h,4 + h,oh,2h,b) 

/i 3 2 - (2/1,1/1,3 + hflh,A)h,2 + /i,oA 2 3 + ^iJiA 

_l_ ( — /l,4^1,2 + ^1,3-^1,2 + h,lh,5h,2 + hfill,4, ~ ^1,1-^1,3^1,4 _ /l,o/l,3/l,5) Z 
/f,2 — 2/1,1/1,3/1^ — /l,o/l,4/l,2 + /l,o/l,3 + /l,l/l,4 

/f,3 ~~ 2/1,2/1,4/1,3 — /l, l/l,5/l,3 + h, 1/1,4 + /l,2^1,5 
/f,2 ~~ 2/1,1/1,3/1,2 — /l,o/l,4/l,2 + /l,o/l,3 + /l,l/l.4 

To determine the orthogonal polynomials for q > 2 we use the property 
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where P^(rii,z) are the orthogonal polynomials for q = 1 with respect to the weight Wi(x). For 
example for q = 2 and j = 2, Eq. (i55l) yields 



P(i,M) 

P(2,M) 



^2 _|_ (h,2h,o — h,oh,2)z h,2h,i — h,ih,2 
h,oh,i — h,ih,o h,ih,o — h,oh,i ' 



z 3 + 



"-^2,3-^1,1 + ^1,4-^2,0^1,1 + ^1,3^2,1-^1,1 — ^1,2-^1,3^2,0 _ -f 1,0^1,4-^2,1 + hfih, 2-^2,3) 



+ - 



+ - 



-^2,2-^1 — -^1,3-^2,0-^1,1 + -^1,2-^2,0 + h, oh, 3-^2,1 ~ A, 2 (-^1,1-^2,1 + A,0^2,2) 

-^2,0-^1,3 ~~ -^1,1^2,2-^1,3 — A,0-^2,3A,3 — A, 2-^1,4-^2,0 + -^1,0-^1,4^2,2 + A, 1-^1,2-^2,3) 2 
-^2, 2-^1 — ^1,3-^2,0-^1,1 — h,2h,lh,l + -^1,2^2,0 + ^1,0-^1, 3-^2,1 ~ ^1,0-^1,2-^2,2 

^2,3-^1,2 ~~ -^l,4^2,lA,2 + -^1,3-^2,1 + -^1,1-^1,4-^2,2 ~ A,3 (A, 2-^2,2 + A, 1-^2,3) 



"-^2,2^1 + -^1,3-^2,0-^1,1 + h,2h,lh,l — -^1,2^2,0 _ ^1,0^1,3^2,1 + h,oh,2h,2 

4.2 Lax operators 

Let us introduce dressing operators Gj according to 

/ 0i (s, z) = (G; &)(s, z), G; := ^ a in (s) If, a in (s) := (a n+1 ) 0i (s), 



n>0 



where s 6 T// and £j(s, z) := z Si . In the matrix representation they are given by the triangular 
matrices 



G, 



( Goo G01 G02 \ 

G11 G12 G13 . . . 
G22 G23 



1 G n m — C^i,m— n(^) 



Si=—m 



\ ■ ■ ■ J 

The corresponding inverse operators are characterized by expansions of the form 

1 1 



n>0 



aio{s) a i(s) 



We define the Lax operators Z{ by 



Zi := Gi T* GJ 1 . 



It follows at once that they can be expanded as 



n>0 



where 



7i(s) = a i0 (s) (T*(3 i0 )(s) = Vi(s). 



(58) 



(59) 



(60) 
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Proposition 8. The functions /oj satisfy the equations 

zf 0i (s,z) = (Zif i)(s,z) + a 0i (s)5 Si o, Vs E T H . 
Proof. From the definition of Gi we have 

zf i(s,z) =Gi(z£i) = Gi (T* (&) + S Si0 ) = {Zif 0i )(s,z) + a i0 (s) 5 Sz0 , 
where we have taken into account that 

T t n (5 SiQ ) = 5 s< - n ,o = °> Vn > 1, a € 17/. 



(61) 



□ 



4.3 The second system of string equations 

The diagonal solutions 

*(z) =^<j) a {z)E a 

a, 

of the condition (j^J) corresponding to the function g(z) of (}4"6"]) are characterized by 

<9 2 w f + 0o u>i - 4>i Wi = 0, i = 1, . . . , q. 

Hence, setting cfio = we get 

$(z)=Y,V'(c i ,z)E i . 

i 

The corresponding covariant derivative is 

D z f:=^-J2V(c i ,z)fE i , 



so that we may write 



D z ( E f) = djoo E + J2( d ^i ~ Via, z) fa) E ( 

i 

In order to take advantage of the last identity we observe that (jBTj) can be generalized to 



n-l 



z n fa(s, z) = (Z? f oi )(s, z)-J2 Pfi, I) ( s > <Wr.o, Vs E T n . 

r=0 

» 



(62) 
(63) 
(64) 

(65) 



where the coefficients r \(s, z) are polynomials in z. On the other hand we have that as z — > oo 



((T*) n f 0a )(s,z) = { 



C{ — ) ;■-<>. for n - U. £ -/?. 



c(-) z Si , for a = i / j and Sj < — n, > n — ^' 



(ir/ oi )(s,z) = 



0, for Sj > — n 



z s \ ra>0. 



(66) 



?n+l 



We are now ready to prove the following result. 
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Proposition 9. The function f satisfies the equation 



(D z + H) (E f)(s, z) = J2 Ai(s, z) E 0l , Vs e T n , 



i=i 

where 7i is the operator 

i=i 

and Aj(s, z) are functions of the form 

Ai(s,z) =^2p( itn) (s, z) 8 Si+ni0 . 

n=l 

with P(i tn )(s, z) being polynomials in z and iVj = degreeV(ci, z) — 2. 
Proof. From (|64p it follows that 



D z (E f)(s,z)+H(E f)(s,z) 



00 



-En 



^ 2 /0i + ^ V '( C 3> Z j)(j,+) foi ~ V '(°i: z ) hi 
3=1 



Eoi- 



o, -0(-J **\ 



Using ([UfJ]) we find 

d z foo + ^2V(c j ,Z j ) ijr 

(y"(ci,Zi) ( i i+) - v'(ci,zj) fa = (y'i^Zi) - v'(ci,zj) fa + e>(~^ 

On the other hand (|65p implies 

n>l 

where 

n-2 

Ai(a,2f) := ^nc in ^p^ 1 ' (s, z) S Si+r>0 . 

n>l r=0 

Hence Eq. ([70|) says that 

« 1 
(£> z + W) (£7 f)(s, z)-J2 Ai(s, *) ^ = O(-) /o(«, *)■ 

Z J *y 



1=1 



The first member of this equation is a solution of the Riemann-Hilbert problem for all 
that from Proposition 2 the statement follows. 



As a consequence we deduce the string equations 
Theorem 4. The multiple orthogonal polynomials of type II verify 

Q 

d z P(n, z) + ^2 V'{cj, %)(,-.+) P ( n > *) = °- ( 74 ) 
3=1 

4.4 Orlov operators 

We define the Orlov operators .Mi by 

M t :=Gi-{s t -l)-T t -GT l . (75) 
They satisfy [2j,A^j] = I and can be expanded as 

A4i = X>n(*)Zf, (76) 

n>l 

where 

Uj(s - ej) 

Proposition 10. The functions fa satisfy the equations 

d z fa(s, z) = (M t fa)(s, z), Vs G F n . (78) 
Proof. From the definition of G{ we have 

d z fa = Gi (( Sl - 1) z" 1 &) = G • («i - l)CZi &) = Gj • ( a< — 1) • Tj • GrVw = M /«• 

□ 

5 The large-n limit 

The large-n limit of multiple orthogonal polynomials is closely connected to the quasiclassical limit 
of the functions fo a (s,z). In this section we will consider these functions for large values of the 
discrete parameters s» 

Si >> 1, Vi (Type I case); s» << — l,Vi (Type II case). 



Note that in particular the string equations (153j) and (|67j) simplify since all the 5 terms vanish. As a 
consequence the resulting equations are the same for both types of multiple orthogonal polynomials 
and can be summarized as follows: 

z fa = [Ti + Ui(s) + 52 j Vj(s) T*^j fa, Va, i; 

(79) 

k d z fa = -H fa, d z fa= + V'ici.ZM fa. 
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In order to define the large-n limit we introduce an small parameter e, define slow variables 

Q 

ti'—esi,; t := - } ]tj, t := (h, . . . ,t q ), (80) 
i=i 

and rescale the exponents of the weight functions (llip and (|46f) as 

Wi{e, z) = exp y =F J , 

where the exponent sign is negative (positive) for polynomials of type I (type II). Moreover, we 
perform a continuum limit in which as e — > 0, the discrete parameters Sj tend to +oo (— oo) for 
the type I case (type II case) and t a become continuous variables. 

The problem now is to determine solutions fo a (e,t,z) of (|79p defined for t on some domain 0, 
of W 1 , that have the quasiclassical form [16] 

f 0a {e,t,z) =z So <*- 1 exp^S a ), S a = t a log z + ^ — S^, (81) 

6 n>0 Z 

where 

S Q n = e fe S<$ (t), n > 0; § o = 0. 

fc>0 

Note the leading behaviour 

f 0a (e,t,z) =z & ^~ 1 exp(is Q + 0(l)), as e - 0, (82) 

where 

S a (t,z):=t a log z + J2^ S *n(t), S an =S®, S 00 = 0, (83) 

n>0 

are the classical action functions. 

In terms of slow variables the operators Tj and T* become translation operators 

T i = exp(-ed i ), T* =T~ 1 =exp(edi), % := (84) 

Hence, we have the following useful relations 

If 1 /oa = exp ( =f diS a + 0{e)) f 0a . (85) 

It is now a simple matter to deal with the corresponding dressing and Lax-Orlov operators. Indeed, 
expressing the functions (|8ip in the form 

h„ = (*>„ + Y, ^-jr^) »p (7 l °s -*) • 

n>l 
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we have 



foi = Gi exp ( — log z J , Gi :=^2 a in (e, t) Tf 



n>l 



Zi-. = GiTr l G;\ Mi-Gi-ti-Ti- Gr 1 
We can also introduce Lax-Orlov associated with foo- In fact we may do it in q different ways 

foo = G$ exp log z) , Gf = 1 + J] a „(e, t) IT", 



n>l 



(0 ._ rrW 



In terms of Lax-Orlov operators and taking into account the assumption (|81jl the system of 
string equations ([79]) becomes 



(86) 



z/o a = ^ a /o a = m + u,-(e,t) + ^ v fc (e,t)T£) /oc, Va,j; 



e d z foo = Mo foo = -H foo, e 9 2 / oj = Mj / 0j - = ( - H + V( Cj , ZM /«, 



for all choices Z = Z$,Mn = MX'. It follows from (ED that the recurrence coefficients u, and 



f (0 



l o 



Vj can be written as quasiclassical expansions of the form 



u ? - = u 



(87) 



n=l 



n=l 



5.1 Leading behaviour and hodograph equations 



Our next aim is to characterize the leading behaviour of the solutions fo a of (|86p . More concretely 
we are going to see how the leading terms 

u := (ui(t),...,u q (t)), v := (vx(t),...,v q (t)), 

of the recurrence coefficients (j87|) are determined by a system of hodograph type equations. 

In order to formulate the classical limits (z a ,m a ) of the Lax-Orlov operators (Z a ,M a ) we 
observe that as a consequence of the first group of string equations in ([86]) we have that 



Ti + \ii(e,t)) f 0a = [Tj + Uj(e,t)) f 0a , Vi,j,a. (88) 

Then, using ([85]) we get 

exp ( - di S a (t, z) ) + Ui(t) = exp ( - dj S a (t, z) ) + Uj(t), (89) 
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In view of these identities we define z a (t,p) by the implicit equations 

p = exp ( - di S a (t, z a (t,p))) + Ui (t). (90) 

Notice that according to (|59"|) these definitions are independent of the value of the index i used in 
(1901) . Moreover, ([90]) imply 

d i S a (t,z a ) = -log(p-u i (t)). (91) 



From the asymptotic expansion (|83p of the action functions S a and the defining equations (|90p it 
is straightforward to prove that the Lax functions can be expanded as 



Vo n (t) 



p^cc, 



n=l 



(92) 



Vi{t) 
p - Ui(t] 



+ ^2vin(t) (p ~ Ui(t)) n , p 



n=0 



On the other hand, we define the corresponding Orlov functions m a (t, z a ) by 

m a (t, z a ) := d z S a (t, z a ). 



(93) 



The definitions (|90p and ()93[) provide the classical limits of the Lax-Orlov operators. Indeed, from 
(i86j) it follows at once that 



(Z a f 0a )(t,z a (t,p)) = z a (t,p) f 0a (t,z a (t,p)), 



(M a foa)(t,z a {t,p)) = [m a + 0(e)) / 0a (t, z a (t,p)) 



(94) 



for all choices of Zq = Zq\Mo = -Mq' . In particular this means that all the pairs of Lax-Orlov 
operators (Zq\ «Mq ) have the same classical limit given by (zo(t,p),mo(t,p)). 

Theorem 5. The Lax-Orlov functions satisfy the classical string equations 



z = zi 



E(u,v,p), 



m = m 1 - V'(ci,zi) 



where 



E:=p + Y, 



Vk{t) 



m q - V'{c q ,z q ) = -H(u,v,p), 



k=l 



(95) 



(96) 



and ( )(fc,+) stand for the projections of power series in (p — Uk) n , (n £ Z) on the subspaces 
generated by (p — Ufc)~ n (n > 1). 

Proof. Taking into account that 

f 0a (t, z a (p,t)) = ((p - Uj ) n + 0(e)) /oa(t, z a {t,p)), n = ±1, ±2, . . . , 
it is easy to see that the equations ([9"5]) are the classical limit (e — > 0) of the system (|55j) . □ 
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In view of the first group of equations in (|95p , it is clear that the functions u and v are the only 
unknowns for determining the Lax-Orlov functions. However, the Lax-Orlov functions must verify 
the correct asymptotic expansions. Obviously, the functions z a = E satisfy ([92]) . Nevertheless, 
Eq. (|83p requires the Orlov functions to satisfy 



ta ti S an {t) , . 

m a = 2^ — ^+T~ i as z a -> °°, (97) 

Za n>l Za 

and this behaviour must be compatible with the second group of equations in (|95p 

mo = -H(u,v,p), mi = V'(ci,E) - H(u,v,p), (98) 



where we have already inserted the substitutions Z{ = E. Let us analyze the equations ([98]) in 
terms of series expansions as p — > oo for mo, and as p — > Uj(t) for m,;. If we take into account that 

-=0((p-«i)), -^— = 0(1); V'( Ci ,E)-H = 0(1), j^i, p-t^i), 

Zi p- Uj 



then the consistency between (|98j) and (|97j) requires 

jf ^-H(u,v,p) = -t Q , (99) 

and 

' / dp ^(ci.^it.i;,^)) - H(u,v,p) _ Q 



7l 2Z7T p - Uj 

(100) 

dp V'(ci,E(u,v,p)) - H(u,v,p) ^ 



2m (p-m) 2 Vi' 

These conditions are obtained by comparing coefficients of p _1 and (p — Ui(t)) n with (n = 0, 1) in 
the equations ([98j) for tuq and mj, respectively. Here ji are positively oriented small circles around 
p = Ui such that p = Uj is outside 7, for all j 7^ i, and 70 is a large positively oriented circles which 
encircles all the 7$ (see fig.l). 

Identifying the coefficients of the remaining powers p~ n and (p — Ui(t)) n in (|98p determines the 
Orlov functions in terms of (u,v). 

The equation ([99]) is a consequence of the second group of equations in (jlOOp and the fact that 
to := — ti- To see this, notice that 



<z> H(p) dp = j> H(p) dp E(p) dp, 

J 70 J 70 



/ (F'( Ci , £(p)) - F(p)) ^ dp = - I (V'(a, E{p)) - H{p)) d p E(p) dp. 



Hence 



23 



Figure 1: 



/ H( P ) dp + T i (Via, E( P )) - H{p)) Vl dp 

- I H{p) d p E{p) dp-^f (V(ci, E{p)) - H(p)) d p E(p) dp 

J-/0 j J~ii 



I H(p) d p E(p) dp-^i V'(ci, E(p)) d p E(p) dp = 



where we have taken into account that V'(ci, E{p)) d p E(p) = d p ( V(ci, E{p))j . Moreover, H(p) d p E{p) 
is a rational function of p with poles at the points pi = Ui only and 



7o 



^7i~0 in C\{pi,...,p q }. 



Therefore we are finally lead to the system (|100|) of 2q equations for determining the 2q functions 
Ui,Vi. These equations are of hodograph type as they depend linearly on the parameters t and c,. 
For example the first few terms are 



en + 2 Cj 2 Uj + 



[Ci2 ~ C j2 ) Vj 



Ui - Uj 



2 a 



i2 



3^1 



(Ui - Uj) 2 



+ 



+ 

u 

Vi 







(101) 



5.2 Connection with the Whitham hierarchy 

If we assume that the coefficients Cj of exponents of the weight functions (jlip and (|46p are free 
parameters and write them in the form 



{tali 



j "arti ■ ■ ■) 



(102) 
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then, as we are going to see, the solution of (|95p turns out to determine a solution of the Whitham 
hierarchy of dispersionless integrable system [11] . 
Let us introduce the modified Orlov functions 



m a = V'(t a ,z a ) +m a . 
It is clear that (z a ,fh a ) solve the system 



(103) 



Zq = Z\ 



{m = mi 



-q: 



m, 



(104) 



Moreover, they are rational functions of p with poles at the points pi = Ui only. Furthermore, they 
satisfy the asymptotic properties (i92|) and 



m 



nt an z a + - n+l 

n>l n>\ 



as Zr\ 



DC. 



Thus, the functions (z a ,m a ) satisfy all the conditions of Theorem 1 of [15] and, as a consequence, 
they verify the equations of the Whitham hierarchy 



dz 
dt 



fin 



(105) 



fm 



where the Poisson bracket is given by 



IPG} — ——-—— x — t 
' dp dx dx dp ' 

and the Hamiltonian functions are 



n 



/in 



(106) 



log(p-Ui), n = 0, (j, = i = l, 



Here (-)(o,+) stands for the projector on {p n }^ =0 ■ 

In this way we conclude that (z a ,fh a ), as functions of the coupling constants Cj = to — U, 
determine a reduced solution of the Whitham hierarchy. This property is in complete agreement 
with the results of recents works |16j which prove that the universal Whitham hierarchy can be 
obtained as a particular dispersionless limit of the multi-component KP hierarchy. Moreover, as it 
as been observed in [24J, appropriate deformations of the Riemann-Hilbert problems for multiple 
orthogonal polynomials determine solutions of the multi-component KP hierarchy. In fact these 
deformations correspond to the flows induced by changes in the parameters t a . Indeed, for both 
types of multiple orthogonal polynomials, (|102j) implies 



d tan g=[z n E a ,g]. 



Therefore the covariant derivatives 



Dan f ■ 9 an f + Z f E a , 
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are symmetries of the corresponding Riemman-Hilbert problems. Hence, using Proposition 2 one 
concludes that 

d an f + (z n fE a r 1 )^f = 0, (107) 

where ( )_ stands for the projections of power series in z k , (k S Z) on the subspaces generated 
by z~ k (k > 1). The equations (|1U7|) constitute the linear system of the multi-component KP 
hierarchy. 



6 Applications: random matrix models and 
non-intersecting Brownian motions 

As we have seen the multiple orthogonal polynomials of type I are the elements /oi of the funda- 
mental solution of their associated RH problem. Thus, in the quasiclassical limit we have 



Ai(n,z) ~ -exp (^Si(t,z))j, as e -> 0, 



where Si = Si(t,z) are the classical action functions defined in (|83p . Hence 



ed z log Ai(n,z) ~ d z Si(t,z) - ^ = m»(t, z) - |. (108) 



On the other hand, if we denote by X{ the roots of Ai(n, z) we have 

rii—X 



d z log Ai(n,z) = — 



Z X{ 
1=1 

Thus if we assume that in the large-n limit the roots of Ai are distributed with a continuous density 
Pi = Pi{x) on some compact (possibly disconnected) support Jj C M 

e V — / ^ dx, as e -> 0. (109) 

^— ' t: — t-.- Ij.Z — X 



Z Xj 



from (|108p and (|109p we deduce the important relation 

mi ( z ) = f ^-dx, (110) 

JliZ-X 

where mi, 1% and pi depend on the slow variables t. This means that the Orlov functions mi are 
the Cauchy transforms of the root densities p^. Hence, they determines the distribution of roots in 
the large-n limit according to 

m,i + {x) — rrii-(x) = — 2iirpi(x), x 6 Ij. (HI) 

Moreover, from (|97p we see that 

p l (x)dx = t i . (112) 
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On the other hand, the multiple orthogonal polynomials of type II represent the element /oo of 
their associated RH problem. Therefore, in the quasiclassical limit we have 

P(n, z) ~ exp ^— So(t, z)^j , as e — > 0, 

Thus if we assume that in the large-n limit the roots of P(n, z) tend to be distributed with a 
continuous density po = po(x) on some compact support Iq C M, we deduce 

m ( z ) = [ ^Mdx, (113) 
7 Jo 2 - x 

where tuq, i"o an d Po depend on the slow variables t. Thus the Orlov function mo is the Cauchy 
transform of the density po and therefore 

m,o+(x) — mo-(i) = — 2in po(x), x £ Jo- (114) 

Note also that 

/ po(jc) dx = t - (115) 
Jlo 

The string equations (|95p provide also useful information to determine the limiting supports 
and the root densities. They imply 

m (z) = -H(po(z)), rrii(z) = V'(ci, z) - H{pi(z)), 

where p a (z) denote the q + 1 inverses of the map 

verifying 

Po(z) = z + Pi(z) =Ui + as z -> oo. 

Therefore (jllip and (|114p reduce to 

ff(p a+ (a;)) - iT(p _(x)) = 2in p a (x), x £ I a . (116) 

In general the limiting supports Iq, may consist of several disconnected segments 

d a 

la = lak 
k=l 

which, due to (|116p . constitute the branch cuts of the functions H(p a (z)) . As a consequence the 
end-points of the segments I a k are the branch points of these functions, which are in turn given 
by the critical points x% of the function z(p) = E{p) 

Xl = E( gi ) £ R, d p E( gi ) = 0. (117) 
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6.1 The Hermitian matrix model 

For q = 1 the multiple orthogonal polynomials of type II reduce to the orthogonal polynomials on 
the real line associated to the weight function w = expV(c, z). These polynomials are connected 
to the random matrix model of nxn Hermitian matrices pQ-[2] 

Z n = J dM exp (Tr V(c, AT)) , (118) 

through the crucial relation 

P n (z)=E[det(z-M)}, (119) 

where E denotes the expectation value with respect to the probability measure determined by 
(I118p . This means that in the large- n limit the eigenvales of M converge with unit probability to 
the roots of P n . As a consequence the root density po of the family of polynomials represents the 
eigenvalue density of the matrix model. 

The Hermitian matrix model provides an appropriate example to illustrate all the aspects of 
our method for characterizing the quasiclassical limit. In this case we set e := 1/n, to = 1 and we 
have 

v 

zip) = E(u,v,p) =p-\ . 

p — u 

Here u and v depend on the coupling constants c = (c\, C2, ■ ■ .) and can be determined by means 
of the hodograph equations (f99j) - (|100p 

70 2z7r J 7l 2iir p-u 

By introducing the change of variable p — u — > p these equations are equivalent to the well-known 
system [T] 

/ |P V'(c, p + u + -) = -l, l^- V'(c, p + u + -) = 0, (120) 

which characterizes the spherical limit in the Hermitian matrix model of 2D gravity. Here 7 is a 
large positively oriented circle around the origin. 

The critical points of E are q± = u ± y^u, so that the support of eigenvalues is 

I=[x-,x+], x±:=m±2 v / w. (121) 

We use (|116p to determine the density of eigenvalues according to 

H(p 0+ (x)) - H(po-(x)) = 2mp (x), x e (122) 

Furthermore, the two inverses of z(p) are 

PoO) := -(z + u+ y/(z- x-)(z - x + ) ), pi(z) :=^(z + u- \J [z - xJ)(z - x + ) ) (123) 



and we have 

H(p (z)) = V'(c,z(p))^ +) 



p=p {z) 
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Now, using the identities 

V 2 , 

= z — p, p = {u + z)p — zu — V, 

p — u 

it is clear that there exist polynomials a k (z) and (3 k {z) satisfying 



(z(p) k ) n =a k (z) + p k (z) P0 (z), (z(p) k ) = a k {z) + f3 k {z) Pl (z). (124) 



P=Po{z) v ' (1,+) 



p=pl(«) 



In particular, taking into account that 



p (z) = z + O (^j , pi(z) = u + o(-) ; as z — > oo, 



from (|124p we deduce 



«^) = - =-(-_===) , (125) 



where ( means the projection of power series in z n , (n £ Z) on the subspace generated by 
z n , (n > 0). Hence it follows that 

H iPo( z )) = ^kc k \ a h -. 1 {z) + I5 k -i{z) Pq{z)^ , 
fc>l 

and therefore we get 

= ^y^^c fe /3 fc _i(x)(po+(^) -po-(x)) = --^-f—= V ^ C ^ =) ^{x-x^)(x + - x), 
which represents the well-known eigenvalue density for the Hermitian model in the one-cut case. 

6.2 Gaussian models with an external source and non-intersecting 
Brownian motions 

For q > 1 the multiple orthogonal polynomials of type II are connected to the Gaussian Hermitian 
matrix model with an external source term AM [6j-[9j, where A is a fixed diagonal nxn real 
matrix. The partition function of this model is given by 

Z n = J dMexp^-TrQ^M 2 - AM^. (126) 

It turns out that if the eigenvalues of A are given by dj, (j = 1, . . . , q) with multiplicities rij, then 
the expectation values 

P(n, z) = E[det(z - M)], n := (m, . . . ,n q ), (127) 
are multiple orthogonal polynomials with respect to the Gaussian weights 



Wj(x) = exp(aj x — — x 2 ) 
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These matrix models are deeply connected to one-dimensional non-intersecting Brownian motion 
[20]-[22]. More concretely, the joint probability density for the eigenvalues (Ai, . . . , A n ) of M is the 
same as the probability density at time t E (0, 1) for the positions (x\, . . . , x n ) of n non-intersecting 
Brownian motions starting at the origin at t = and forming q groups ending at q fixed points 
bi, (i = 1, . . . , q) at t = 1. The corresponding dictionary for this duality is 



Xi = — = , CLu — b. . 

3 v^T^J Vi-t 

We discuss next an example of application to the large-n limit of non-intersecting Brownian 
motions. Let us consider an even number n non-intersecting Brownian motions ending at two 
points ±6 with ri\ = n,2 = n/2 [9]. In this case the slow variables take the values t\ = = —1/2. 
Moreover, we have 



z 2 z 2 I t 

V(ci,z)=az , V(c2,z) = — az , a: 



2' 2' Vl-r 

and 

z( P ) = E{p)=p+^^ + ^—, H{p) = ^_ =p _ z (p). (128) 

P — U\ P — U2 p — U\ P — U2 

Using the hodograph equations (jlOlh one finds 

1 

ui = a, u 2 = -a, vi = v 2 = -, 

so that 

p A — a 1 

The corresponding algebraic function p = p(z) satisfies the Pastur equation [23] 

p 3 - zp 2 + (1 - a 2 )p + a 2 z = 0, 

which defines a three-sheeted Riemann surface. The restrictions of p(z) to the three sheets are the 
functions p a {z) characterized by the asymptotic behaviour 

Po(z) = z + Pi(z) = Ui + o(^J, i = 1,2; as z -> oo. 

There are four critical points of z(p) which give rise to four branch points ±xi, ±X2 in the z-plane 
where 

Vl + 8a 2 + 3 Vl + 8a 2 -3 

X1 = Q1 Vl + 8a 2 + l ' X2 = Q2 Vl + 8a 2 -l ' 



- + a 2 ± \\Jl + %a 2 . 



Qi,2 w 2 2 



It is easy to see that x\ is real for all a > 0, while X2 is real for a > 1 (x2 < x±) and purely 
imaginary for < a < 1. Now, from (|122p and taking into account that H(p) = p — z(p) we deduce 
that the eigenvalue density is given by 

Po(x) = — (H(j> 0+ (x))-H(j> _(x)) = — ( Po+ (x)-po-(x)), xel . (129) 

2 ITT 2 ITT 
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Figure 2: Limit support for Brownian motions with two symmetric endpoints for 6 = 1 

Po Po 




-3-2-10123 -3-2-10123 



Figure 3: The density of Brownian motions po(x) for a = 1/2, 3/4, 1 and 3/2, respectively 



Using Cardano's formula for po one finds 

2x 2 + 6 (a 2 - 1) - \^(r(x) - yV(x) 2 - 4s(x) 3 ) 2 3 

PoOs) = , , 

2 5 /3 y/3 vr 3/ r ( x ) _ v / r (s) 2 -4s(s) 3 

where 

r(x) := -2 x 3 + 18 a 2 x + 9 x, s(x) := x 2 + 3 (a 2 - 1). 
The form of the support Io depends on the analytic properties of the function po(z) ( see [6]- [9]): 

a) For < a < 1 the function po is analytic in C — [— xi,xi] and Iq = [—xi,xi]. 

b) For a > 1 the function po is analytic in C — ([— x±, —X2] U [X2, xi]) and Iq = [ — ^1, — ^2] U [#2, ^l]- 
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